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RADAR  CROSS  SECTIONS  OF  DISTRIBUTED  TARGETS  LOW  OVER  WATER  IN 
THE  PRESENCE  OF  MULTIPATH  REFLECTIONS 


INTRODUCTION 

Recently  there  was  some  interest  in  determining  the  detection  range  of  small,  distributed  targets  (such 
as  periscopes)  low  over  the  ocean  in  the  presence  of  multipath.  Average  radar  cross  sections  (RCS)  of 
those  objects  measured  in  sea  trials  routinely  exceeded  those  measured  in  free  space  by  5  to  7  dB.  In  order 
to  explain  this,  the  RCS  of  a  line  of  equally  spaced  reflectors  suspended  quasiperpendicularly  low  over 
the  ocean  is  treated  in  the  following  calculations. 


Theory 

Figure  1  shows  the  four  reflection  paths  from  a  point  source  over  the  ocean.  The  height  of  the  radar 
above  the  water  is  denoted  as  H,  the  height  of  the  point  reflector  is  x,  the  horizontal  distance  between 
radar  and  target  is  R,  and  the  reflection  coefficient  on  the  water  surface  is  p.  For  horizontal  polarization 
and  sea  state  0  (i.e.,  no  waves),  the  reflection  coefficient  p  =  -1.  For  higher  sea  states,  p  tends  towards 
zero,  more  so  for  shorter  electrical  wavelengths  and  higher  grazing  angles.  Figures  2,  3,  and  4  show  |  p  | 
as  a  function  of  rms  waveheight  with  grazing  angle  as  parameter  for  3,  10,  and  35  GHz,  respectively 
(adapted  from  Miller  et  al.  ). 


The  lengths  of  the  return  paths  are  (see  Fig.  1) 


•  Direct  path 

•  First  reflected  path 

•  Second  reflected  path 

•  Third  reflected  path 


2  si 

s  1  +  s3  +  52 
s2  +  s3  +  si 
2(s2+s3) 


The  voltage  received  by  the  radar  is  the  vector  sum  v  of  the  voltages  from  all  four  paths  (the  path  loss 
due  to  the  fourth  power  of  range  was  taken  out): 


v  =  <*/*.  si  +  2p  e'  2lt/x(sl+s2+s3)  +  p2  e>4,t/x(s2+s3))  v/" 


(1) 


where  v/~is  the  voltage  that  would  be  created  from  the  target  in  free  space. 

si  =  -J(x  -  H)2  +  R2 


(2) 


A.R.  Miller,  R.M.  Brown,  and  E.  Vegh,  “New  Derivation  for  the  Rough-Surface  Reflection  Coefficient  and  for 
the  Distribution  of  Sea-Wave  Elevations,”  IEE  Proceedings  131(H2),  April  1984. 


Manuscript  approved  August  20,  2001. 
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Fig  3 . —  Specular  reflection  coefficient  over  sea  at  10  GHz 


For  (x-H)1  «  R\ 


si  *  R  +  1/2  (x-H)2  /R 


Letting 


s2  +  5*3  ~  R  +  1/2  ( x+H)  /R 


s  1+52+53  =  2 R  +  (x2+H2)/R 


a  =  2n/X  (x+H)2/R  b  =  2n/X  (x-H)2/R  c  =  2n/X  (x2  +  H2)/R, 
v  =  gWwx  +  2p  dc  +  p2  g/n)  v/ 

Iv|  =  |(/>  +2p^  +  pV")l  | v/| 


The  square  of  the  absolute  value  of  v  is  obtained  by  multiplying  v  by  its  conjugate  complex.  It  can  be 
shown  that 

I  V  | 2  =  (  1  +p2  +2pcos(»))2  I  v/l2,  where  u  =  4n  3  H  /  (XR).  (6) 
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Fig.  4  —  Specular  reflection  coefficient  over  sea  at  35  GHz 


Determination  of  Free  Space  RCS  of  a  Point  Target  from  the  Multipath  Maxima  and  Minima 

If  the  multipath  minima  and  maxima  of  an  encountered  point  target  are  known  as  it  moves  at  constant 
height  above  the  water  surface  towards  or  away  from  the  radar,  the  following  calculations  are  used  to 
determine  its  free  space  RCS. 

The  maximum  occurs  when  p<0  and  u  =  {2k+\)n,  where  k is  an  integer.  Then,  with  Eq.  (6), 

I  v  1 2max  =  /w  =  (  1  +  I  P  I  )4  I  v/| 2  (7) 

I  v  1 2min  =  pmin  =  (1  - 1  p  I  )4  I  vf\ 2  (8) 

Let  w  =  |  v  |  max  /  I  v  Lin  the  ratio  of  the  maximum  to  minimum  voltage  observed  by  the  radar  receiver 
as  the  point  target  moves  from  a  multipath  maximum  to  the  subsequent  minimum  at  constant  height,  then 
the  free  space  radar  cross  section  of  the  target  Gf  can  be  determined  by  the  formula 

<Jf=amax(  1  +  2/Vw  +  1/w)  / 16 ,  (9) 


where 

tfmax  =  RCS  of  point  target  measured  at  a  multipath  maximum 
w  “  !  ^min 

Gmin  =  RCS  measured  at  adjacent  multipath  minimum 
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The  cross  sections  must  be  corrected  to  be  referenced  to  the  same  range  by  using  the  1/R4  law  of  the  radar 
equation. 

Multipath  Pattern  from  a  Single  Reflector 

Figure  5  shows  the  familiar  multipath  pattern  resulting  from  a  point  target  moving  to  or  from  the 
radar  at  constant  height  of  1  m  above  the  sea  level.  The  free  space  return  is  referenced  to  0  dB.  For  clarity, 
the  R'4  law  of  the  radar  equation  was  omitted  in  the  result,  representing  the  normalized  RCS.  The  peaks 
are  12  dB  above  free  space.  The  reflection  coefficient  p  was  assumed  to  be  -1,  i.e.,  sea  state  0.  The 
frequency  was  10  GHz,  the  height  of  the  radar  above  sea  level  was  59  m. 


Fig.  5  —  RCS  with  multipath,  single  reflector  at  1  m  above  sea  surface  (radar  height:  59  m;  sea  state:  0 
(calm);  p  =  -1 ;  and  f  =  1 0  GHz) 


Multipath  Pattern  of  Distributed  Targets 

The  target  is  assumed  to  be  flat  rising  from  the  water  surface  to  a  height  xmax.  It  is  tilted  from  the 
perpendicular  by  an  angle  a.  According  to  Huygens  principle,  the  target  can  be  represented  by  many 
equally  spaced  reflectors,  provided  that  the  spacing  is  small  in  comparison  to  the  electrical  wavelength 
(see  Fig.  6).  Then,  using  Eq.  (5),  we  can  compute  the  normalized  radar  return  from  this  target  by  simply 
summing  the  voltages  resulting  from  each  reflection  point: 

nmax 

|  v |  =  |z  (^b  +  2p  e!c  +  p2 ef*  )  I  \vf\  /nmax,  (10) 


with 
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b  =  2n(xv  +H) 2  /(X  R )  +  4k  xv  /X  sin(a)  (10a) 

The  added  term  in  b  is  taking  into  account  the  additional  phase  shift  due  to  the  tilt.  The  same  term  is 
added  to  c  and  a. 


H 


R 


alpha  =  tilt  angle  of 

reflector  array 
from  vertical 


txmax 


* 


* 


Fig.  6  -  Geometry  of  the  reflector  array 


Equation  (10)  was  used  to  obtain  the  following  results.  Figure  7  shows  the  pattern  resulting  from 
three  reflectors  spaced  perpendicularly  above  the  sea  surface  at  0,  0.5,  and  1-m  heights.  Some  of  the 
minima  are  filled  in,  and  the  maxima  have  come  down  from  12  dB  to  about  6  dB. 

Figure  8  shows  the  same  scenario  as  Fig.  1  but  with  667  reflectors  spaced  equally  within  zero  to  1  m 
height  perpendicularly  above  sea  level.  The  result  may  be  somewhat  unexpected.  The  multipath  pattern 
has  disappeared;  the  return  is  independent  of  range  equal  to  about  6  dB  above  free  space  peak  return  of  a 
linear  array. 


Figure  9  shows  the  scenario  of  Fig.  8  but  with  no  multipath,  i.e.,  free  space.  This  represents  the 
reflection  pattern  of  a  linear  array,  with  the  aspect  angle  of  the  array  varying  with  the  range.  At  a  range  of 
1000  m,  the  RCS  of  the  particular  sidelobe  peak  is  23  dB  below  the  free  space  maximum  and  29  dB 
below  the  return  seen  with  multipath  as  in  Fig.  8.  It  is  clear  that  multipath  plays  a  decisive  role  in  this 
case.  The  sea  surface  and  the  reflector  act  like  a  dihedral  reflector. 

Next,  we  consider  what  happens  when  the  reflectors  are  tilted  off  perpendicular  by  3.3  degrees,  with 
no  multipath  reflections  present.  At  a  range  of  1000  m  we  see  the  peak  return  of  the  flat  plate,  because  at 
that  range  and  59  m  height  of  the  radar,  the  beam  hits  the  reflector  array  perpendicularly.  Figure  10  shows 
the  result.  The  maxima  and  minima  are  due  to  the  sidelobes  of  the  reflection  pattern.  When  we  add 
multipath,  we  see  the  result  shown  in  Fig.  1 1.  The  peak  at  1000  m  range  is  almost  unaffected  by  multipath 
reflections,  but  at  ranges  below  500  m,  the  return  is  increased  by  up  to  20  dB. 
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RANGE  (m) 

Fig.  7  • —  RCS  with  multipath,  three  reflectors  at  0,  0.5,  and  1  m  perpendicularly  above  sea 
surface  (radar  height:  59  m;  sea  state:  0  (calm);  p  =  -1) 
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Fig.  8  —  RCS  with  multipath,  667  reflectors  perpendicularly  above  sea  surface  from  0  to  1  m 
(xmax  =  1  m;  xmin  =  0  m;  H  =  59  m;  tilt  angle  a  =  0;  sea  state:  0  (calm);  p  =  -1;  and  f  =  10 
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Fig.  11  —  RCS  with  multipath,  667  reflectors  above  sea  surface,  tilted  3.38°  from 
perpendicular,  from  0  to  1  m  (xmax  =lm;  xmin  =  0  m;  H  =  59  m;  tilt  angle  a  =  3.38°;  sea 
state:  0  (calm);  p  =  -1;  and  f  —  10  GHz) 


The  angle  of  3.38  degrees  was  interesting  because  it  rendered  the  perpendicular  aspect  angle  at  1000- 
m  range.  However,  larger  tilt  angles  such  as  15  degrees  may  be  of  more  interest. 

Figure  12  shows  the  return  of  the  reflector  array  tilted  15  degrees  off  perpendicular  away  from  the 
radar,  with  multipath  present.  Again,  at  about  220  m  range,  the  peak  return  of  the  flat  plate  pattern  is  seen, 
because  at  that  point,  according  to  the  geometry,  the  radar  beam  hits  the  array  perpendicularly. 
Remarkably,  there  is  no  enhancement  of  this  peak  due  to  multipath,  because  for  that  geometry,  the  main 
lobe  of  the  reflection  pattern  of  a  flat  plate  does  not  hit  the  sea  surface.  With  increasing  height  of  the 
radar,  this  peak  moves  out  farther  in  range. 

Figure  13  shows  the  same  RCS  as  in  Fig.  12,  but  with  no  multipath.  It  is  seen  that  at  very  close  range 
like  at  100  m,  the  RCS  is  enhanced  by  multipath  by  about  8  dB,  at  900  m  by  about  5  dB.  But,  what 
happens  at  ranges  farther  out? 

Figure  14  is  like  Fig.  13,  but  for  a  range  out  to  8000  m.  The  pattern  seen  is  not  due  to  multipath  but 
rather  to  the  varying  aspect  angle  of  the  array.  Figure  15  is  like  Fig.  14,  but  with  no  multipath  reflections. 
It  is  shown  that  the  peaks  are  increased  by  5  to  6  dB  due  to  multipath  reflections.  Figure  16  shows  the 
pattern  for  a  reflection  coefficient  of  -10  dB  (0.32).  Hardly  any  difference  is  seen  between  no  multipath 
and  a  -10  dB  reflection  coefficient. 


RCS  (dB) 
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RANGE  (m) 


Fig.  12  —  RCS  with  multipath,  667  reflectors  above  sea  surface,  tilted  15°  from  perpendicular, 
from  0  to  1  m  (xmax  =  1  m;  xmin  =  0  m;  H  =  59  m;  tilt  angle  a  -  15°;  sea  state:  0  (calm); 
p  =  -1;  and  f  =  10  GHz) 


Fig.  13  —  RCS  with  no  multipath,  667  reflectors  above  sea  surface,  tilted  15°  from 
perpendicular,  from  0  to  1  m  (xmax  —  1  m;  xmin  =  0  m;  H  =  59  m;  tilt  angle  a  =  15°;  free  space, 
no  multipath;  and  f  =  10  GHz) 
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RANGE  (m) 

Fig.  16 —  RCS  with  multipath,  667  reflectors  above  sea  surface,  tilted  15°  from  perpendicular, 
from  0  to  1  m  (xmax  =  1  m;  xmin  =  0  m;  H  -  59  m;  tilt  angle  a  =  15°;  with  little  multipath; 
p  =  -0.32;  and  f  =  10  GHz) 


Distributed  Reflector  with  Rough  Surface 

The  influence  of  the  roughness  of  the  reflector  is  considered  next.  The  reflector  is  made  up  of  an  array 
of  small  reflectors  located  along  a  line  from  the  sea  surface  up  to  a  height  of  1  m,  the  line  being  tilted  off 
the  vertical  by  the  angle  a.  Now,  every  adjacent  10  reflectors  are  moved  out  of  the  line  towards  the  radar 
by  delta_r,  delta_r  being  equally  distributed  between  zero  and  A,/4.  Every  10  adjacent  reflectors  had  the 
same  delta_r.  There  were  20  reflectors  per  wavelength. 

Figures  17  and  18  show  the  effect  of  reflector  roughness  for  a  smooth  sea  (p  =  -1)  and  for  free  space 
(p  =  0),  respectively.  Figures  19  and  20  show  the  same  for  the  range  extended  to  8  km.  It  is  remarkable 
that  the  peak  at  the  range  where  the  radar  beam  hits  the  reflectors  perpendicularly  is  very  little  degraded 
from  its  original  value  of  0  dB.  For  ranges  greater  than  1000  m,  the  RCS  with  full  multipath  is  greater  by 
5  to  10  dB  than  the  free  space  RCS.  The  maxima  and  the  minima  resulting  from  the  reflection  pattern  of 
the  smooth  array  are  averaged  out  for  the  rough  array  at  ranges  above  1000  m,  in  that  geometry. 

CONCLUSION 

If  the  surface  of  the  reflector  is  rough  (i.e.,  the  rough  surface  peaks  are  smaller  than  X/4),  then  the 
maxima  and  minima  resulting  from  the  sidelobes  of  the  array  reflection  pattern  are  smoothed  out  as  a 
function  of  range.  The  RCS  of  the  array  in  full  multipath  condition  is  5  to  10  dB  larger  than  that  of  free 
space  for  ranges  of  1000  m  or  greater,  for  the  geometry  considered. 
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Fig.  19  —  RCS  with  multipath,  667  reflectors  above  sea  surface,  simulating  a  rough  reflector 
tilted  15°  from  vertical,  height  from  0  to  1  m;  H  =  59  m;  with  multipath;  p  =  -1;  and  f  =  10  GHz 


0  1000  2000  3000  4000  5000  6000  7000  8000 

RANGE  (m) 

Fig.  20  —  RCS  without  multipath,  667  reflectors  above  sea  surface,  simulating  a  rough  reflector, 
tilted  15°  from  vertical,  height  from  0  to  1  m;  H  =  59  m;  without  multipath;  p  =  0;  free  space;  and 
f=  10  GHz 
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With  a  reflection  coefficient  of -0.32  (-10  dB),  there  is  practically  no  difference  between  free  space 
and  the  multipath  environment  for  this  scenario.  With  smooth  sea,  sea  state  0,  reflection  coefficient  of  0 
dB,  the  peaks  in  the  l-to-7-km  ranges  are  enhanced  by  5  to  6  dB.  The  peaks  are  not  due  to  multipath  but 
rather  to  the  reflection  pattern  of  a  linear  array  as  a  function  of  aspect  angle.  As  the  range  changes,  so 
does  the  aspect  angle. 


Appendix 

MATLAB  PROGRAMS  FOR  FIGURES  2  THROUGH  4  (SPECRHO.M) 

AND  5  TO  16  (DISMULT.M) 
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%  Program  name  specrho.m 

%  Calculates  the  specular  reflection  coefficient  rho 

%  as  a  function  of  wave  height,  frequency  and  grazing  angle  as  parameter 
%  According  to  Miller,  Brown  and  Vegh, 

%  IEE  Proceedings  Vol  131  Pt.  H  No.  2  April  1984 
%  f  =  frequency  (GHz) 

%  sigh  =  rms  wave  height  (m) 

%  alf  =  grazing  angle  (degrees) 

%  rho  =  absolute  value  of  specular  reflection  coefficient  (dB) 


f =3 ; 

lam  =  .  3  /  f ; 
alf = [0 . 1  0.2  0.5  1]  ; 
alfarc=alf *pi/180 ; 
p=sin(alfarc) / lam; 
delsig=0 . 001 ; 
sigmin=0 . 1 ; 
sigmax=10 ; 
m=l ; 

while  m<5 
n=l  ; 

sig=sigmin; 

while  sig<sigmax 
sig-sig+delsig; 
g=sig*p (m)  ; 
u  ^  2* (2*pi*g) ^2 ; 
if  u>300 
u=3  00 ; 
end 

y=  20* ( loglO (bessel j ( 0 , j  *u) ) -u*logl0 (exp ( 1 ) ) ) ; 
if  y<-30 
y  =  -  3  0  ; 
end 

rhodb(mfn)=  y; 
sigh (n) =sig; 
n=n+l ; 

end 

n=ly 

sig=sigmin; 

m=m+l 

Grid 

semilogx (sigh, rhodb ( 1 , : ) , sigh, rhodb (2 , : ) , sigh, rhodb ( 3 , : ) , sigh,  rhodb (4 , : ) ) 
grid 
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%  Radar  response  of  distributed  reflectors  low  over  the  sea  surface 
%  with  multipath  reflections.  The  grazing  angle  is  assumed  to  be  small 
%  {smaller  than  the  Brewster  angle)  .  The  Reflector  is  an  assembly  of 
%  vertically  equally  spaced  equal  reflectors,  emulating  a  flat  plate. 

%  The  whole  assembly  is  tilted  from  the  perpendicular  by  an  angle 
%  alfa  of  less  than  20  degrees. 

%  The  reflection  coefficient  rho  on  the  sea  surface  is  equal  to  -1  for 
%  glassy  smooth  sea  (sea  state  0)  for  small  grazing  angles  and  horizontal 
%  polarization.  The  fourth  power  of  range  of  the  radar  range  equation  is 
%  omitted  in  the  result.  The  total  RCS  of  the  reflectors  in  free  space 
%  is  normalized  to  0  dBsm. 

%  Parameters : 

%  lam 
%  h 
%  rho 
%  xmax 
% 

%  alf 
%  nrefl 
%  r 
%  rmin 
%  rmax 
lam=0 . 03 ; 
h=59  ; 

%xmax=l ; 
xmax=l ; 
xmin=0  ; 

%alf=2 . 78 ; 

%alf =3 . 38 ; 
alf =15 ; 

%alf =0 ; 

nrefl=  20* (xmax-xmin) /lam; 

%  The  spacing  of  the  reflectors  is  lambda/20 
%nref 1=2 ; 

%x=  0 : xmax / nre  f 1 : xmax ; 
x=xmin: (xmax-xmin) /nrefl : xmax; 

%x=l : 1 : 1 ; 

%x=xmax : xmax : xmax ; 

%  rmax=8000; 

%  rmin=200; 
rmax= 1000; 
rmin=l; 

%rmin=670; 

%x=0 . 3 ; 

rho=- 1 ; 

%  rho=0? 

%rho=-0 .  316  ; 
nmax=1000 ; 
n=l  ; 

alfarc=alf /180*pi; 
talf =tan (alfarc)  ; 
while  n<nmax 

r=rmin+  (n-1 )  /nmax*  ( rmax- rmin)  ; 

ar(n)=r; 
k=2*pi/lam; 
phi=2*k*x*talf ; 
aa=k/r* (x+h) . *2  +  phi; 
bb=k/r* (x-h) . ^2  +  phi ; 
cc=k/r* (x. ~2+h~2 )  +  phi ; 

v  (n)  =abs  (sum(exp(  j  *bb)  +2* rho* exp  ( j  *cc)  +rho/v2*exp(  j*aa)  )  )  / length  (x)  ; 
if  v(n)  <  le-3 
v(n) =le-3 ; 
end 

lv(n) =20  * loglO (v(n) )  ; 


=  Electrical  wave  length  (m) 

=  Height  of  radar  above  the  water  surface 
=  reflection  coefficient  of  the  water  surface 
=  maximum  height  of  distributed  reflectors  (m) 

(the  lowest  reflector  is  assumed  to  be  on  the  water  surface) 
=  tilt  angle  of  the  reflector  assembly  from  perpendicular 
=  number  of  reflectors  (depends  on  lam  and  xmax) 

=  range  (m) 

=  minimum  range 
=  maximum  range 


